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1 Introduction 



The hadronic width of the Z° boson, I\ = T(Z° — > hadrons), or the ratio 

T h ^ G F Ml 

which is closely related to the familiar ratio R in e + e~ annihilations into hadrons, provide 
theoretically very clear conditions for the verification of perturbative QCD. This is due 
to several favourable circumstances: 

• nonperturbative power corrections are expected to be negligible at the scale Mz 

• perturbative corrections are calculated up to next-to-next-to-leading order (NNLO) 

• the dependence of perturbative expansion for (JJ) on the choice of the renormalization 
scheme (RS) is weak 

• in calculating QCD corrections to (|]) the running a s can be taken as that corre- 
sponding to five effectively massless quarks, af> 

• the effects of the finite bottom quark mass on the expansion coefficients of I\ in 
this couplant have been calculated up to the NLO 

• the explicit dependence of I\ on the top quark mass m t has been calculated in the 
large m t expansion up to the NNLO 

• the statistics of the data is large, and thus the experimental accuracy rather high, 
in particular compared with the closely related process of e + e~ annihilations into 
hadrons. 

On the other hand, at the scale Mz the magnitude of QCD corrections to the basic 
electroweak decay mechanism is smaller than for quantities at lower energy scales, such 
as the r-lepton semileptonic decay width and therefore more difficult to pin down. In 
this part of the paper we give quantitiative estimates of some of the above-mentioned 
uncertainties and address the related question: assuming the validity of QCD and taking 
into account these uncertainties, how accurately can the basic QCD parameter a s be 
determined? 

As any meaningful discussion of the quantitative importance of higher order QCD 
corrections depends on the ability of experiments to 'see' them, we start by recalling the 
relevant experimental data Jl||: 

M z = 91.187 ± 0.007 GeV (2) 
T = 0.99528 ± 0.00023 GeV (3) 
T h = 1.7407 ±0.0059 GeV, (4) 

which imply the following relative errors, relevant for further discussion 

Ar 4 Ar h AR Z , 

— - = 2.3 10" 4 , — - = — - = 3.5 10~ 3 . (5) 

1 1 h 
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Note that the precision Ar of the determination of T is more that an order of magnitude 
better that that of Th and can thus be neglected with respect to AIV 

Because the dominant part of perturbative QCD predictions for Th has the same 
generic form as does Rz, i.e., 

(6) 

and a s (M z )/ii = 0.037, these errors allow a s (M z ) to be determined to within about 
8.5% acccuracy. Translated into the sensitivity to A, this amounts to a factor of 1.9 
uncertainty. Improving further the accuracy of the data by a factor of two would allow it 
to be extracted with an error of only 36%. 

2 The RS dependence: general considerations 

Over the last 15 years the problem of the renormatization scheme dependence of finite 
order approximants to perturbation expansions in QCD (and other theories) has been the 
subject of lively and sometimes even heated debate. From time to time a 'resolution' of 
this problem is announced, but invariably it turns out that these 'solutions' contain the 
original ambiguity in some guise ot another. We intend to provide a concise and balanced 
review of all the various approaches to this problem, but emphasize at the very beginning 
that, in our view, there is no clear winner. Nevertheless, as the dependence of finite order 
perturbation expansions on the choice of RS is a very real phenomenological problem, 
which cannot be ignored, we think the right question in this context is: How sensitive are 
these approximations to the choice of RS? But even in this question there hides a catch, 
as to give it a concrete meaning we first have to define the set of 'allowed' RS. The point 
is that without some restriction on the considered RS we could get essentially any result 
we want. But again, as the selection of the 'allowed' RS is inevitably a subjective matter 
and may, moreover, depend on the quantity in question, the best we can do is choose a 
couple of approaches which are sufficiently general, have some rationale behind them and 
define the theoretical 'error' with respect to this set of RS. 

Whether the theoretical error of some quantity should be considered large or small is, 
of course, not given a priori, but depends on the accuracy of experimental data to which 
it is compared. These experimental errors for the quantities related to the Z° decay were 
estimated in the preceeding section. 

In this section only the RS dependence of physical quantities will be discussed. For 
unphysical quantities, such as the Green functions with anomalous dimensions, the situ- 
ation is more complicated and some of the approaches are not directly applicable. This, 
however, is not a serious limitation, as what we are actually interested in are clearly only 
the physical quantities. 

Furthermore, we shall consider only the case of QCD with rif massless quark flavours. 
The reasons for this restrictions are twofold. First, the relations resulting from the renor- 
malization group (RG) considerations and expressing the internal consistency of the renor- 
malized perturbation theory, have yet to be worked out for the general massive case. The 
lack of such relations precludes the general quantitative discussion of the RS dependence 
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problem, possible in the massless case. In the Z° mass range, however, we can with great 
accuracy consider QCD with five effectively massless quark flavours. This statement is 
quantified in Section 3. 

In the following subsections we first discuss the quantitative description of the freedom 
connected with the choice of the RS ('kinematics' of the RS dependence problem) and 
then briefly review several of the approaches to choosesing one of these RS ('dynamics' of 
the RS problem). We emphasize this distinction, as the two aspects are frequently mixed 
up. Only the latter aspect is really of substance, the former being merely a matter of 
convention and bookkeeping. 

2.1 The description of the RS dependence 

Consider the generic perturbation expansion for the physical quantity of the form 

r(Q) =a(RS) [r + rx(Q, RS)a(RS) + r 2 (Q, RS)a 2 (RS) + ■ ■ •] ; r„ = 1, (7) 

which appears in the expression for Th and ([!]). Q in (]?]) denotes generically some external 
momentum on which r depends [] and a(RS) is the renormalized couplant a = a s /ir (the 
adjective renormalized will be dropped in the following). There are many different ways to 
quantify the dependence of such physical quantities on the RS. As a matter of convention, 
we shall adopt the one suggested in Ref. 0. First, we should define the meaning of the 
renormalization scheme itself. In massless QCD there must be some parameter with the 
dimension of the mass, for the moment loosely denoted as A, that sets the basic scale 
of the theory. Once this parameter is given, any quantity can in principle be calculated 
as a concrete number. For a given A, fixing the RS means specifying the values of all 
perturbative coefficients rfc(RS), as well as the value of the expansion parameter a(RS) 
itself. 

The labelling of the RS suggested in Ref. starts with the familiar equation 
da(/i, RS) 



(3(a) = —ba 2 (fi, RS) (l + ca(/i, RS) + C2d 2 (/i, RS) + 



dln/i 

expressing the dependence of a on the scale /i, which inevitably appears in the theory 
during the process of renormalization. The first two coefficients on the r.h.s. of (|]), i.e. 
b, c are unique functions of the number rif of massless quarks 

, llN c -2n f 51N c -19n f , n . 

b = -; c = -, (9) 

6 ' IIN C -An f ' V ; 

but all the higher order ones are completely arbitrary. Once they are given and some 
initial condition on a is specified, (|8|) can be solved. The way of specifying the boundary 
condition is ambiguous, but its choice is a matter of convention only. One way of doing 
this is via the scale parameter A introduced in the following implicit equation for the 
solution of (181) 101: 



, i Z 1 1 i ca f j 
Mn-^ = — hem h / ax 

A a 1 + ca Jo 



1 1 

+ 



x 2 B^)(x) x 2 {l + cx) 



(10) 



1 In view of the application to Th we restrict our considerations to quantities depending on a single 
external momentum Q. 
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where 

B {n \ X ) EE (1 + CX + C 2 X 2 + ■■■ + Cn^X* 1 - 1 ). (11) 

Note that as the integral on the r.h.s. of ( JTD ) behaves like 0(a s ) the higher order coef- 
ficients Ck]k > 2 have no influence on the value of A. Note also that the parameter A 
introduced in ([10 ) differs from the A used in most phenomenological analyses by a factor 



close to unity: A = A (2c/(V) C//b . At the next-to-leading order (NLO) - i.e., keeping only 
the first two terms in ( pTOj ) - the solution of ( ]TD| ) is often approximated by the first two 
terms of its expansion in powers of inverse logarithms ln(/i/A): 

/ ,n 1 c ln(ln(/i 2 /A 2 )) 

^ /A) = MnWA)-^ InVA) + - (12) 

The dependence of the couplant a on the parameters % > 2 is determined by equa- 
tions similar to (||) [0: 

da(/i, Cj) _ fa*+ 2 

=A = - /5(a) io (^ d "' (13) 

which are uniquely determined by the basic /3-function in (§) and thus introduce no 
additional ambiguity. 

It is obvious that a unique definition of a(fi) at some /J requires, as well as the spec- 
ification of the coefficients c^i > 2, the specification of the boundary condition, - i.e., 
for instance, the value of A. It is convenient to introduce the concept of the renor- 
malization convention (RC), which is associated with a fully defined solution of (§): 
RCee {A,q; i > 2}. As, however, \i always enters this solution in the ratio ///A, we 
can either: 

• select one of the solutions to (§), which we call referential renormalization convention 
(RRC) and vary // only or 

• fix fi by identifying it with some external momentum - for instance Q - and vary 
the solution of @, i.e., for fixed coefficients q the value of A instead. 

Both these options are completely equivalent and it is merely a matter of taste as to which 
one to use. We prefer the former. To vary both simultaneously is legal, but obviously 
redundant. 

As the choice of the RRC is a matter of convention only, we cannot associate any 
physical meaning to the scale \x itself. It serves to label the RS, but only in a given RRC. 
In two different RRCs the same \i may correspond to different values of the couplant 
a as well as the coefficients r^. We emphasize this point as in many papers the RS is 
chosen by identifying /j with some 'natural' physical scale of the process, such as the 
external momentum Q. Although such a natural scale can usually be identified, its mere 
existence does not help fix the arbitrary scale /x, as to get a unique RS, the RCC also has 
to be specified. This is usually tacitly assumed to be the MS, but there is no theoretical 
argument for this choice, except that in this RRC the coefficients are often small. If, 
however, the magnitude of the coefficients of the perturbative series for physical quantities 
would be the criterion, we would be naturally drawn to the effective charges approach, 
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described below, where they actually vanish. In other words, because the choice of the 
MS as the RRC is merely a convention, there is no reason to set \x = Q. 

The above relation ( |TUD allows the expresssion of ///A in terms of a and q and thus 
the labelling the RS by means of the set of parameters a, q; i > 2. Using this way of 
labelling the RS is very convenient as there is then no need to introduce the RRC and 
also no possibility of referring to the 'natural' scale to fix the RS. 

In the NNLO order - i.e., taking into account also the first nonunique coefficient c 2 - 
we have the equation 



bint 
A 



i ca 

- + c In == 

a v 1 + ca + C2d z 



(14) 



where 

/(a,ca) 



2c 2 - c 2 



d 



2c 2 



d 



2c 2 <2 + C 
1 

2c 2 a + c — d 



arctan ■ 



In 



arctan - 
a 



c — d 



d 



4c 9 - c 2 ; 4c? > c 2 



-In 



c + d 



d = Jc 2 - 4c 2 ; 4c 2 < c 2 .(15) 



2c 2 a + c + d 

Its solution depends on the value of c 2 . We distinguish three different cases: 

• c 2 = (resp. Ci = 0,i > 2), defining the so called 't Hooft RC 

• c 2 > 0, when /3(a) < is monotonously decreasing function of a and the situation 
is therefore qualitatively the same as for c 2 = 

• c 2 < 0, when (3{a) has the infrared fixed point at a*(c 2 ), given by the equation 
/3[a*(c 2 )] = 0. The corresponding solution of (H) then approaches /iimte value at 
\i = and consequently 



limr (3) (Q) 



s (c 2 ) (l + n(/i = Q)a*(c 2 ) + r 2 (c 2 ,/i = Q)a* 2 (c 2 ) + • • ■) (16) 



has a finite infrared limit at the NNLO. This case is discussed in detail in Refs. 

§!• 

Note that the possibility of an infrared stable limit of finite order approximants, which 
starts at the NNLO, does not have to survive the incorporation of still higher order 
corrections and its physical relevance is therefore questionable. 



While the explicit dependence of the couplant on q is given in ([14]) and (|i~5|), the 
dependence of the coefficients on them is determined by the requirements of internal 
consistency of the perturbation theory. They imply that any finite order approximant 



rM(Q)=5> fc , 



N-l 



T(li,c h pui <N -I) 



(17) 



k=0 



must satisfy the following consistency conditions: 

drW _ N+1) drW 



d ln/i 



0(a 



0{a 



(If 
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Iterating these equations we find: 



r^Q/p) = b\n^ + ri ^ = Q) = b\nt-p(Q/~A) 



r 2 {Q/n, c 2 ) = P2 - c 2 + r x + o- x , (19) 

and similarly for still higher orders. In the above relations the quantities p,p 2 etc., are 
RG invariants - i.e., contrary to the coefficients rj~, they are independent of the choice 
of the RS. Note that all the dependence of the perturbative approximants on Q comes 
exclusively through the invariant p(Q/A), which can be written as 

p = Mn(Q/A RRC ) -n(jji = Q,RRC), (20) 

where the apparent dependence on the chosen RRC actually cancels between the two 



terms in (|20|) 



A nontrivial part of any perturbative calculation boils down to the evaluation of these 
invariants, the rest being essentially a straightforward exploitation of the RG considera- 
tions based on (|l8|) . Substituting for the term 61n(/x/A) in r\ the expression (|l0l) , using 
(|19|) and inserting the resulting into (|17), any finite order approximant r^ N \Q) can be 



expressed as an explicit function of the parameters specifying the RS, i.e., a, q; i > 2, and 
the invariants pf. 

r (N \Q) = f(pjj < N - 1; a, a, i < N - 1). (21) 

In this representation the RS dependence of NLO and NNLO approximants is quanti- 
tatively described by one- and two-dimensional manifolds, respectively, and the problem 
of choosing the RS is equivalent to selecting one particular point on these manifolds. 
Considered as a geometrical exercise we identify certain special points on these manifolds, 
corresponding to stationary points, where the variation of the approximants with respect 
to the free parameters vanishes locally. 

In the next subsection we shall briefly describe some of the criteria for choosing the RS, 
which will define the set of RS for which we shall later estimate the theoretical uncertainty 
of perturbative calculations of the quantities of interest. We shall discuss in some more 
detail the approaches described in subsections 2.5 and 2.6, as there have recently been 
some new developments in them. 



2.2 Fixed RS calculations 

Because of the computational simplicity and explicit gauge invarince, all the multiloop 
calculations are nowadays done using the dimensional regularization technique. Within 
this technique the MS renormalization prescript ionP] is often preferred on the grounds 
that it absorbs in the definition of the renormalized couplant the terms proportional to 
In 4-7T— je, which are considered to be artefacts of the dimensional regularization technique. 
In our way of labelling the RS, MS corresponds to definite values of all the coefficients 
rk, Ck and a fixed, but numerically undetermined, value of a, which must be extracted 

2 By 'prescription' we mean the specification of the scale fi (by identifying it with some natural scale 
Q) as well as of finite parts of all counterterms necessary to cancel the UV divergencies. Specifying the 
prescription implies the specification of the RS in the above-defined sense, but not vice versa. 
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from comparison with experimental data. This choice of the RS is very commonly used in 
phenomenological analyses, but there is no obvious reason why it should be preferred to, 
for instance, the MOM-like RS or any of the choices discussed in the following subsections. 
In geometrical terms this is reflected in the fact that the corresponding point on the 
hypersurfaces defined in (^) occupies no special position. 

2.3 Principle of Minimal sensitivity (PMS) 

In this approach, suggested in Ref . |§ , the RS is fixed by demanding that 



da dc,- 



0, (22) 



i.e., the N-th order partial sum has locally the property that the full expansion must 
satisfy globally. Though there is in general no quarantee that such a stationary point is 
unique or exists at all, in practical applications to lowest order QCD quantities it works. 
At the NLO, when only a labels the RS, (p2|) reduces to 



ca f ca \ , 

2 - 2pa + 2caln + ca = 0, 23 

1 + ca Vl + ca/ v ' 

and its solution has the form opms = (Vp)[1 + C( ca PMs)]- 

At the NNLO we have two coupled equations for derivatives of A 2 \a, c 2 ) with respect 
to a and c 2 , which must be solved numerically. For the quantity (^) such a stationary 
point exists for any p if p 2 < and for p > p m i n (p2) if P2 > ||. 

Note that, contrary to the fixed RS approach, the PMS selects the RS which depends 
on the type and kinematics of the process under study. The same holds for the methods 
discussed in the next three subsections. 



2.4 The method of effective charges (ECH) 

The basic idea of this approach || is to choose the RS is such a way that the relation 
between the physical quantity and the couplant is the simplest possible one. For the 
quantity (0) it means: 

r{Q) = a ECH . (24) 

In this approach there is no problem with the convergence of the perturbation expansion 
© itself, but it reappears in the perturbation expansion of the corresponding /^-function 
(see below). 

The conditions under which the parameters a, Ci, or p, Ci, can be chosen in such a way 
that (^) holds can be read directly off the consistency conditions flT8|). At the NLO, ( f2~3| ) 
implies the following equation for a ECH 

1 1 CCL 

— hcoln = p, (25) 

a 1 + ca 

which has a solution Oech = (l/p)(l + C(caECH)), differing from a PMS merely by the term 
of the order C(coech)- At the NLO the values of oech thus correspond to intersections 
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of the curves defined as 



r (2) (a,c 2 ) =a[2- pa + caln — V (26) 
V f + ca ) 

with the straight line r^ 2 > = a. For p > there is always just one such intersection in the 
physically relevant range a > 0, while for p < there is none in this range. 

At the NNLO the situation is more complicated, as the condition ([24]) does not by 
itself determine uniquely oech, but merely implies the equation 

ri(a) + r 2 (a, c 2 )a = 0, (27) 

which has, depending on the value of c 2 , either no solution, or one or two solutions, giving 
different giech- Going to still higher orders the ambiguity of this approach grows even 
further as new free parameters, q, crop up. In || this ambiguity is avoided by demanding 
that each of the coefficients 7*j vanishes individually. Assuming this restricted version of 
the ECH method - i.e., demanding = 0, i > 1 - we get the following expression for the 
associated /3-function: 

^ flECH = (3 E CH(a) = -bal CH \l + pia EC H + P2a E CH H 1 > ( 2 §) 

amp, L J 

where pi = c. The coefficients cech,i of the ECH /3-function /?ech thus concide with 
the RG invariants pi, introduced in fll9"D. To express oech as a function of the external 
momentum Q, it is convenient to write it as a solution of the equation 

1 . . ca EC H Q {on , 
+ cln— = 61n 7 , ( 29 ) 

^ECH 1 + CaECH ^ECH 

where Aech defines the 'effective' A parameter, associated with the quantity under study. 
It is related to Ars in any fixed RS simply as 

A E CH = A RS exp(r 1 ( Ai = g,RS)/6). (30) 

As the ECH approach seems to offer a very simple and natural 'solution' to the RS 
problem one might naturally ask: where has all the ambiguity discussed in subsection 2.1 
actually gone? In fact it has not disappeared entirely and reemerges, as discussed in Ref. 
H, in a somewhat disguised form, even within the ECH approach. 



2.5 The method of Brodsky, Lepage and MacKenzie 

This method || borrows its basic idea from QED, where the renormalized electric charge 
is fully given by the vacuum polarization due to charged fermion-antifermion pairs. In 
QCD the authors of this method suggest fixing the scale \x with the requirement that all 
the effects of quark pairs be absorbed in the definition of the renormalized couplant itself, 
leaving nothing in the expansion coefficients. In the case of the quantity (0) and up to 
the NNLO, 



r(Q) = a(fi, RRC) 1 




71/T-ii 




a(fi, RRC)^J , (31) 
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where we have now written the rif dependence of the coefficients explicitly, it amounts 
to the requirement that \i be chosen in such a way that rn(fi/Q, RRC) = 0. The problem 
with this 'scale-setting' method is that the resulting scale as well as r BLM depend on 
the choice of the RRC! The is due to the fact that for a given [i the separation of the 
coefficient r% into the two parts A + rifB is not unique, but depends on the RRC used. 
This problem exists in principle QED as well, but there the quark loop effects can be 
rather unambiguously absorbed in the renormalized electric charge via the MOM RRC, 
which for massless quarks gives the same B as the MS-like ones. This is no longer true 
in QCD, where various types of MOM-like RRC in general give different values of B, 
different again from that of MS-like one 0. 

As emphasized in the general discussion above, fixing the scale without also simulta- 
neously fixing the RRC does not, however, determine the RS, because the choice of the 
RRC is equally important as that of the scale fi. The resulting ambiguity of the BLM 
approach, pointed out a long time ago ||, has not yet been satisfactorily resolved. In 
practical applications one usually starts from the MS RRC. 

It is claimed in Ref. || that the BLM-improved expressions for the physical quantities 



have small NLO coefficients. However, as shown in Ref. [10], this is not necessarily the 



case when the BLM approach is generalized to higher orders. Let us first recall the main 



steps of the generalization suggested in Ref. [TIJ 



Within the class of the co called 'regular' RC the n/ dependence of the expansion as 
well as /^-function coefficients is polynomial in n/: 

ri = r w + r u n f 

r 2 = r 20 + r 21 n f + r 22 n 2 f 

(3 = b = boo + boin f 

ft = bc = p 1 o + Pi 1 n f (32) 

j3 2 = bc 2 = (3 20 + (3 21 n f + f3 22 n 2 + p 23 n 3 f 

(3 2 = bp 2 = (3 20 + (3 21 n f + (3 22 n 2 } + j3 23 n 3 f . 
Note that the scheme-invariant coefficient f3 2 contains the n 3 term as observed in Ref. 



The generalization of the BLM approach suggested in [10] assumes that the chosen 



scale /i is determined by the following perturbative expansion: 

V 2 = /4LMU +7l( n /) a 0"BLM) + (33) 

where /iIlm i s given in and 71 = 710 + 7n^/- The parameters 710 and 711 are process 
dependent and can be determined from the following system of equations 3 

/3 23 -/?23 = -2/3 2 l7 ii 



P22 - P22 = -2/9m 7io - 4/3 o/3ii7ii 



3 Note the factor-of-two difference between our definition (g) of the /3-function and that used in Ref. 
PI. 
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hi - An = An fa - rfj - Pur* - 4A,oA)i7io - 2/3 2 o7i 2 i 
A20 - #20 = A)o ( r* 2 - rf J - ^i r* - 2/3o 7i , 



which follows from the general expression for the scheme invariant /^i 

= fa + - /3in - fcr 2 . 



(34) 



(35) 



In the above equations, rjj 1 and r\ are the n/-independent coefficients in the generalized 
BLM procedure. We have already mentioned that in practice the BLM approach is 
applied to the initial series with the coefficients defined in the MS-scheme. Therefore, it 
is necessary to put /5 2 3 = 0. Now consider as an example perturbation expansion for the 
familiar quantity: 



R(s) 



crie 1 e 



hadrons) 



a e T e" 



1=1 



1 + E r k a 

k=0 



k+l 



Applying the above generalization to (p6|) we get |T0 

/4lm = /4isexp(0.69), 701 = 0.11, 711 « 3, 



which implies 
where 



3 E<3 2 



1 + a, + 0.08a: - 23.3a; 



i=l 



i=l 



J2Qi) l-24af 



a* = a /^bliJ 1 + 7i("/)o(^blm) 



(36) 



(37) 



(3* 



(39) 



Notice that the coefficient of the NNLO correction is, indeed, not small. Therefore it is 
not true, as conjectured in Ref. ||, that the BLM-improved perturbative series have in 
general significantly smaller coefficients than the expansions in the MS. 



2.6 RS invariant perturbation theory 

The basic problem of RS dependence can be traced back to the fact that the expan- 
sion parameter, the renormalized couplant a, is not a physical quantity, but rather an 
intermediate variable, allowing us to correlate different physical quantities. As such, it is 
inevitably ambiguous. This problem can be circumvented - at least in part - by express- 
ing one physical quantity directly as power expansion in terms of the other. Consider, 
for instance, two physical quantities admitting the following perturbation expansions in 
some RSf] 

R {1) = a(RS)(l + r! 1} (RS)a(RS) + ---) 

R {2) = a(RS)(l + rf } (RS)a(RS) + ---) (40) 
4 There is no problem generalizing this analysis to the case of different powers of the leading terms. 
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Expressing o(RS) from the first in terms of and substituting into the second equation 
we get 

rW=rU(i + AWrU + ...), (41) 



where A (2>1) = rf } (RS) - r? } (RS), as well as all the other coefficients A^\j > 2 of 
this expansion, are unique. This kind of expansion has already been discussed within 
the so-called 'scheme invariant perturbation theory' in Ref. [L2] and recently resurected 
within the so called 'commensurate scale relations' in |y|. The essence, however, remains 
the same. 

As a special and interesting example of such a relation, consider the one between the 
derivative of ([7p with respect to the external momentum Q and the quantity r(Q) itself. 
It is straightforward exercise to show that this relation reads: 

= -b[r(Q)] 2 (l + cr(Q) + p 2 [r(Q)} 2 • • •) = / ( 42 ) 

where the r.h.s. is nothing else than the 'effective' /3-function introduced above, evaluated 
at r(Q)\ 

In (]4T|) as well as in (42) there is no trace of any RS ambiguity. It is nice to be able to 
show explicitly, as is done in Refs. [12-14], that all the coefficients in these expansions are, 
indeed, RS invariants, but it cannot be otherwise, as they relate two physical variables and 
there is no way their eventual dependence on the RS could be cancelled. Perturbation 
theory serves here merely as an intermediate, but vital, tool for evaluating coefficients 
such as A( 2,J ) and pi. 



If relations such as (f4T|) or (|42| ) are truly unique, do they not solve the whole RS- 
dependence problem? The answer depends on what we expect from the perturbation 
theory. If we are interested merely in relating pairs of physical quantities admitting 
purely perturbative expansions the answer is positive. If, however, we ask the question: 
what is, on the basis of analyses of available experimental data, the QCD prediction for, 
say, the Z° width, then the answer is definitely negative. The point is that in predicting 
the value of I\ from the relations between this quantity and some other physical quantity 
R, we find that the resulting predictions depend on the choice of the quantity R\ What 
has been thrown out the door in the form of the RS ambiguity, comes back through the 
window as the 'initial condition' ambiguity Moreover, this new one is even more 

difficult to handle than the original one! 

We therefore believe that the better way of incorporating the results of QCD analyses 
of many different physical quantities measured in different kinematical ranges is to in- 
troduce some intermediate, no-nunique and thus unphysical, variable, which can then be 
used for QCD predictions of other physical quantities. The renormalized couplant serves 
just this purpose. 



3 Quark mass thresholds in the running a s 

The proper treatment of quark mass effects in the QCD running coupling constant a s 
has so far not been in the forefront of interest of theorists and phenomenologists. This 
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is now changing. Although the quark mass effects are not large, steady improvement 
in the precision of experimental data, coming in particular from a new generation of 
experiments at CERN and Fermilab, combined with significant progress in higher order 
QCD calculations, has led to a renewed interest in quantitative aspects of these effects 



Rig , |17j , |18[ . There are two basic reasons for this, both of which are relevant to the subject 
of this article. 

The first concerns the exploitation of the complete NNLO QCD calculations that have 



recently become available for quantities such as (|36| ) |fl9| or I\ (E2(| |21[ and which exist 
basically for massless quarks only. The NNLO corrections are tiny effects and to include 
them makes sense only if they are large compared to errors resulting from the approximate 
treatment of the quark mass thresholds in massless QCD. To quantify the importance of 
the NNLO correction to the couplant a s , consider the difference between the values of 
a s (MS, Mz) in the NLO and NNLO approximations, assuming five massless quarks and 
taking, as an example, = 0.2 GeV. At the NLO the couplant a = a s /rc is given as a 
solution to Equation (0) for c 2 = and we get a NLO (MS, M z ) = 0.03742. At the NNLO 
and in MS RS, where c 2 (MS,n f = 5) = 1.475, we find a NNLO (MS,M z ) = 0.03665. The 
relative difference between these two approximations, 



Q NLO _ Q NNLO 
^NLO 



= 0.02, (43) 



thus amounts to about 2%. The NNLO correction should therefore be included only if the 
neglected effects can be expected to be smaller than this number. However, as we shall 
see, quark mass threshold effects can in some circumstances be of just this magnitude!^ 
The second reason is related to the problem of comparing the values of a s , determined 
from different quantities characterized by vastly different momentum scales. As recently 
emphasized in an extensive review of a s determinations [p2]| , there is a small, but non- 
negligible discrepancy between the value of a s (MS, M z ) obtained by extrapolation from 
some of the low energy quantities, and ct s (MS, Mz) determined directly at the scale Mz 
at LEP - the latter giving the value higher by about 5-10%. Simultaneously, it has been 
noted in ref. [^] that there is an exception to this behaviour in the case of the ratio 
R T , which, when extrapolated from m T to Mz, gives values of a s close to those measured 
directly at LEP (0.120 ± 0.005 ]2"3fl). The physical relevance of the discrepancy between 



the low energy extrapolations and direct measurements of a s at LEP has very recently 
been emphasized in Ref. In particular, its author argues that the extrapolation 



of a s (m T ) to a s (Mz) is unreliable due to limited control of the power corrections. The 
question of estimating the theoretical uncertainty in the extraction of a s (m T ) from data 
on R T is also discussed in ref. p3 As shown in [18], a part of this overestimate of the 



extrapolated value of a s (M z ) pointed out in 0] may in fact be due to the approximate 
treatment of the c and b quark thresholds. 

We shall now analyze the quantitative consequences of the exact treatment of quark 
mass thresholds at the LO and formulate the conventional matching procedure p6f for 



massless quarks in such a way that its results are so close to those which are exact that 

5 In general, as the magnitude of higher order corrections to a s depends on the renormalization scheme 
employed, so does also the estimate (^). However, if defined as the relative difference between the NLO 
and NNLO approximations, this dependence is not strong (see the concluding paragraph of Section 4). 
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the available NNLO calculation can be consistently included. We shall describe in detail 
the approximation in which the u, d and s quarks are considered massless while the c, b 
and t quarks remain massive. 

As complete multiloop calculations with massive quarks are very complicated and 
available only at the leading order, all higher order phenomenological analyses use the 
calculations with a fixed effective number n/ of massless quarks, depending on the char- 
acteristic scale of the quantity. For relating two regions of different effective numbers of 



massless quarks the approximate matching procedure developed in Ref. jo]] is commonly 
used. It should be emphasized that this procedure concerns only those mass effects that 
can be absorbed in the renormalized couplant. At higher orders there are, however, mass 
effects that remain in the expansion coefficients even after the effects of heavy quarks 
have been absorbed in a suitably defined running couplant. 

In QCD with massive quarks the renormalization group equation for the couplant 
a(fj,) formally looks the same as in massless QCD. The only, but important, difference 
concerns the two lowest order /3-function coefficients, b, (3\ = be, which are no longer 
unique as in massless QCD, but may depend on the scale /i. While in the class of MS-like 
renormalization conventions b = 11/2 — n//3, exactly as in massless QCD, in MOM-like 
ones it becomes a nontrivial function of the scale /J ||27|| : 



b(ii/mi) = \- - - ^ h i( x i), x i = ( 44 ) 

where the sum runs over all the quarks considered, rrii are the corresponding renormalized 
quark masses^ and the threshold function h{x) is given as Ref. p7| . The shape of the 



function h(x) is actually not quite unique and depends on the vertex chosen for the defini- 
tion of the renormalized couplant. This fact was first noted in Ref. |28| and subsequently 



explained in Ref. [£{J. The form of h(x) used below corresponds to quark-gluon-quark 
vertex with massless quarks, using any momentum configuration and any invariant de- 



composition in MOM-like RS [p9|| . It is appropriate for most of the extrapolations from 
low energy quantities to LEP energy range. The same form is valid for the ghost-gluon- 
ghost vertex. The three gluon vertex gives somewhat different form of h(x), though its 



behaviour for small and large x is the same. According to [271: 



j / \ . 2 /"I z 2 [l-z) 2 1 6 12 VIT^ + x . x 2 

h(x)=Qx 2 dz \ , ' — - = 1 H In , = -. (45) 

V 1 Jo l + x 2 z{l-z) x 2 x^VaTx 2 VITx^-x 5 + x 2 V ; 

The last, approximate, equality is a very accurate approximation of the exact form of 
h(x) in the whole range x G (0,oo). This allows a simple treatment of the quark mass 
thresholds at the LO. There is, unfortunately, no analogous calculation of the next (3- 
function coefficient, 6c, for massive quarks. This is one of the reasons why most of 
the phenomenological analyses use the so called 'step' approximation, in which at any 
value of one works with a finite effective number of massless quarks, which changes 
discontinuously at some matching points Consequently, b(nf) effectively becomes a 
function of /x, discontinuous at these matching points, as shown in Fig. la. 

6 For the purposes of this discussion quark masses can regarded as constants. In numerical estimates 
they are identified with the running masses at the scale Mz- 
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exact b(yU.) 
/c = 2.24 
58 




3.5 4 

log, (/x/GeV) 



Figure 1: a) 6(/i/m;) together with three step approximations, corresponding to matching 
at the points \x% = nrnf, i = c,b,t with k — 1, 1.58, 2.24; b) the ratio R a for the same three 
approximations, made to coincide at [Aq — 1 GeV. 



The matching points are assumed to be proportional to the masses of the corresponding 
quarks, /ij = Km;. In principle, a different quark threshold could be associated with a 
different k, but for simplicity's sake we take them equal. The free parameter k, allowing 
for the variation of the proportionality factor, turns out to be quite important for the 
accuracy of the step approximation. At the LO the matching procedure then consists of 
the following relations at the matching points fii (the numbers in the superscript define 
the corresponding effective number of massless quarks)[]: 

/ A (3) \ V36(4) 

C(WA (3) ) = C /A (4) ) ^ A (4) = A (3) ^ (46) 



Km c J 

■ A ( 4 ) \ V36(5) 



C(WA (4) ) = CWA (5) )^A( 5 ) = AWr- (47) 



a™* [Kin, 



A (5) 



1/36(6) 



/A< 5 >) = aJ£ 6 («™t/A (6) ) A (6) = A (5) — . (48 



app v '"<■/ ) app 



Note that each of the intervals of fixed rif is associated with a different value of the A- 
parameter, A*W. The resulting dependence of a(fx/mi) on [i is thus continuous at each 
of the matching points, but its derivatives at these points are discontinuous, reflecting 



7 The parameter A appearing in this as well as the other formulae in this section is the leading order A 
parameter, A LO , which cannot be associated with any well-defined RS. As there is no danger of confusion, 
we shall drop the superscript 'LO' for the remainder of this section. 
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the discontinuity of the step approximations to b(/j,/rrii). This procedure can easily be 
extended to any finite order. Let us point out that the more sophisticated procedure 
for matching the couplants corresponding to different effective n/ developed in Ref. []3"m 
coincides in the LO with the above relations (48-50). To estimate the errors in a s resulting 
from the above-defined approximate treatment of quark thresholds we merely need solve 
the LO equation with exact explicit mass dependence as given in 



(49) 



dln/x 

For our purposes the approximation h(x) = x 2 /(5 + x 2 ) is entirely adequate and yields 

a(/i) = , (50) 

11 3\, /i 1 _ , Jn 2 + 5m 2 



2 3 



A(3) 3 j£,t y/(A<$f + 5m? 



where the fraction | comes from the sum over the three massless quarks u, d and s and A^ 
is the corresponding A-parameter appropriate to three massless quarks. For the heavy 
quarks c, b and t we take in the following m c = 1.5 GeV, m& = 5 GeV, m t = 170 GeV. The 
distinction between the 'light' and 'heavy' quarks is given by the relative magnitude of m; 
and A, the latter being defined by the condition 5m? 3> A. For the above values of m c , rrib 



and m t this condition is very well satisfied. Consequently, for \i <C m i} i = c,b,t (|50| ) 
approaches smoothly a LO for nf—3, while for \i ^> and neglecting A^ 3 ^ with respect 
to 5m?, it goes to 

a{n) ~ 1 - 1 (51) 
,,M, At 1 / \Thm c y/5m b y/Em t \ 6(6) ln 

6(6)ln Al3) + 3 ln (^^^J A(6)( ^ } 
where the parameter A^(/c) depends in general on k and 

, , ^ , / A (3) A (3) A (3) \ MleT ( \ \ W) , s 

A( 6 )(v/5) = A( 3 ) = 4= A( 6 )(l) (52) 

V V5m c V5m b V5m t J \VO/ 

coincides with A( 6 ) defined via the subsequent application of the matching relations (48- 
50) for /t = y/E = 2.24. Even though from the point of view of the matching procedure, 
k is not exactly fixed, the value k = a/5 will be shown to be in a certain sense the best 
choice. The relation between and A® depends nontrivially on k. 

In Fig. la, b(fi) is plotted as a function of \i for the above-mentioned masses of c, b and 
t quarks, together with its step approximations and corresponding to three different values 



of k = 1, a/5, J 5/2. There is hardly any sign of the steplike behaviour of the function h(x) 
in the region of the c and b quark thresholds and only a very unpronounced indication 
of the plateau between the b and t quark thresholds. The step approximations are poor 
representations of the exact h(x), primarily due to the rather slow approach of h(x) to 
unity as x — > oo. However, there is a marked difference between the three approximations. 



16 



While the step approximation with the conventional choice k — 1 understimates the true 
h(x) in the whole interval displayed, and would do so even when some smoothing were 
applied, k = a/5 gives clearly much better approximation as the corresponding curve is 
intersected by the exact h(x) at about the middle of each step. 
In Fig. lb the /i dependence of the ratio 

p _ a ex > (aO /c-o\ 

between the above exact solution (|51~1) and the approximate expressions for the above 
mentioned values of k, is plotted assuming A® = 200 MeV0. As we basically want to 
compare the results of different extrapolations starting from the same initial fio, used 
in the approximate solutions has been rescaled by the factor 1.004 with respect to A^ 
in (|50D , thereby guaranteeing that all expressions coincide at fiQ = 1 GeV. Any deviation 
from unity in Fig. lb is then entirely the effect of an approximate treatment of the heavy 
quark thresholds. Figure 1 contains several simple messages. 

The approximate solutions based on the matching procedure defined in (48-50) are 
generally much better immediately below the matching point than above it, and worst at 
about 5m matc h. This reflects the fact that the function h(x) vanishes fast (like x 2 ) at zero 
but approaches unity for x — > oo only very slowly. Moreover, in the Mz range the effect 
of the c quark threshold is essentially the same as that of the b quark and both are much 
more important than that of the top quark, although Mz/m c « 60, Mz/rrib ~ 18, while 
Mz/mt^ 1/2! 

The effect of varying k is quite important, in particular with respect to the c and b 
quark thresholds. In general, k > 1 improves the approximation above, but worsens it 
immediately below the matching point it. The choice k = a/5, suggested by the asymptotic 
behaviour of (|5T|), is clearly superior in practically the whole displayed interval fi G (1, 10 4 ) 
GeV and leads to an excellent (on the level of 0.1%) agreement with the exact solution 
in this interval. On the contrary, the conventional choice k = 1 leads to a much larger 
deviation from the exact result, which exceeds 2% in most of this region. This discrepancy 
is of the same magnitude as the effects of NNLO corrections to the couplant itself. It 
thus turns out that the effect of an exact treatment of the quark mass thresholds for the 
extrapolation of a s from the scales around fi ~ 2 GeV up to fi = Mz is as important as 
that of the NNLO correction to a s and must therefore be taken into account whenever 
the latter is considered and compared with a s determined at these vastly different scales. 

On the other hand, the preceding discussion tells us little about the accuracy of the ap- 
proximation of five massless quarks directly at the scale Mz-, for instance when calculating 
r^. This question will be addressed in section 5. 

The analysis of quark mass effects in a s presented above strictly speaking holds only 
for the LO. Nevertheless, as both the mass effects and the higher order perturbative 
corrections are small effects, it seems reasonable to expect that the conclusions drawn in 
this section will have more general validity. 

8 The resulting ratio R a (Mz) depends only weakly on A^ 3 ). Note, however, that the current analysis 
of the LEP data gives the value of A^ 3 ) in the range 500-700 MeV in disagreement with the QCD sum 
rules analysis (24|. 
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4 Application to 



We now come to the quantitative estimate of the theoretical uncertainties of perturbative 
QCD predictions for T h . In the preceding section we discussed the approximation in 
which a s is given by an expression corresponding to five massless quarks. In this section 
we quantify the uncertainties resulting from the RS ambiguity, discussed in Section 2. 
The nontrivial dependence of I\ on m;, and m t , coming from effects not included in the 
running couplant corresponding to five massless quarks, is discussed in the next section. 
Nevertheless in order to avoid unnecessary repetition we include the dependence on the 
ratio Mz/m t already in the formulae quoted below and obtained recently in Refs. [ |19| , 
The dependence of the expansion coefficients on the bottom quark mass 

We do not write it out here explicitly, as this 

becoming 



20, GIL El 



has extensively been studied in Ref. [33 



would further complicate the structure of fl54]) , each of the coefficients in (p4|) 
a different function of the ratio m^jMz- At the end of this section we shall merely recall 
the leading contributions to Rz coming from m^jMz terms and discuss their numerical 
importance. 

In order to quantify the theoretical uncertainty related to the choice of the RS, we first 
define the set of 'allowed' RS. As emphasized above, this is to large extent a subjective 
matter. Based on our previous experience we define as a measure of this uncertainty the 
difference between the results obtained (for the same Aj^g-) in the three principal methods 
set out in Section 3: PMS, ECH and MS. This choice is to large extent arbitrary, but as 
the MS RS is used in most phenomenological analyses, we adopt it for the lack of anything 
better. The formulae quoted below are taken from Ref. . 

The basic quantity of interest, Rz, defined in ([!]), has a nontrivial structure which 
mixes the effects of electroweak interactions with those of pure QCD. It can be written in 
the following decoupled form (i.e., for five massless flavours and the explicit m t dependence 
of the expansion coefficients) as the sum of three terms with different electroweak factors 
33] and separated further into four possible combinations of vector, axial vector and 



singlet, nonsinglet contributions: 




(54) 

where gv,i = — 2Qi sin 2 9w, 9a,i = ^3 is the third component of the weak isospin 
and the sums over the electroweak coupling constants equal 

5 



r i = {9v,i + 9% 



44 
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+ — sin 9w — ^ sin 9 
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1.6807 ±0.0012 



(55) 



r 2 = I £ 9V,i 



sin 2 9 



0.42850 ±0.00028. 



(56) 
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The last equalities in the above relations correspond to the world average sin 2 9yy = 
0.2329 ± 0.0005. Note that (|55|) implies 

— - = 6.5 10~ 4 < — - = 7.1 10" 4 < — -, (57) 

r2 m ± h 

which, combined with @, means that the theoretical uncertainties of QCD predictions - 
in the square brackets in ([54]) - should be compared to the error of Rz itself. 

The expansion coefficients entering the above formula can be expressed as functions 
of the ratio x = M z /m t (Mz), where m t (Mz) is the renormalized, 'runnning' mass of the 



top quark, taken at the scale Mz- For rif = 5 and in small x expansions we have |21 



n = 1.409 

+ [0.065185 - 0.014815 lnx]x 

+ [-0.0012311 + 0.00039683 In x]x 2 

+ [0.000061327 - 0.000023516 lnx]x 3 + 0(x 4 ) (58) 

r 2 = -12.767 

+ [-0.17374 + 0.21242 lux - 0.037243 In 2 x]x 

+ [-0.0075218 - 0.00058859 lnx + 0.00038305 In 2 x}x 2 

+ [0.00050411 - 0.00012099 lnx + 0.000031419 In 2 x]x 3 + 0(x 4 ) (59) 

s 3 = -0.41318 

4° p = 0.027033x + 0.036355x 2 + 0.00058874x 3 + 0(x 4 ) (60) 
t 2 = -3.0833 + In x + 0.086420x + 0.0058333a; 2 + 0.00062887a; 3 + 0(x 4 ) (61) 
t 3 = 18.654 + 1.7222 lnx + 1.91671n 2 x 

+ [-0.12585 + 0.28646 lnx - 0.011111 In 2 x)x 

+ [-0.0031322 + 0.012117 lnx - 0.0011905 In 2 x]x 2 

+ [-0.00088827 + 0.00047262 lnx - 0.00017637 In 2 x]x 3 + C(x 4 ). (62) 

Because of different electroweak factors in front of them, each of the expressions in the 
square brackets of fl54D is separately from the point of view of QCD RS invariant. As the 
optimization according to either the PMS or ECH methods does not commute with the 
operation of addition, the first question we have to answer is the order of these operations. 
In the absence of uncertainties in the values of the electroweak factors, the proper way 
would be first to sum all three terms in fl54|) and then to fix the RS. In reality, however, the 
errors of electroweak factors induce uncertainties in values of the coefficients multiplying 
powers of the QCD couplant. To optimize, in one way or another, the resulting QCD 
perturbative expansion ( |5"4] ) in such circumstances is not a well-defined exercise and we 
therefore have chosen to follow the opposite route. We believe that to get an estimate 
of the RS dependence of QCD calculations this second route is adequate. In this section 
we shall discuss the numerical importance of the first two terms in (|54"D for m t = 0. The 



effects of non-zero m t} as well as the third term in (54), the existence of which is also 
closely related to non-zero m t , are dealt with in the next section. The second term of (|54"D, 
multiplied by (|56"D , is given at the LO only and no optimization is possible. Fortunately 



it contributes, in MS RS, a mere 10 to Rz and is thus clearly negligible. 



19 



r f =5, p = 27, p 2 = 



-1 5.45 




The dominant contribution to Rz comes from the first term in (|54f) . The term in the 
square brackets can be written as 1 + r(M^), where r{Q) has exactly the form of (0) 



and, moreover, for massless quarks coincides with the above expression for fl3~6|) [19] . For 
rif = 5 the crucial RG invariant p 2 = —15.45, which implies (for detailed discussion see 
Ref. H]) that the saddle point of the NNLO approximant r^ 3 \p) will occur at the point 
(apivis,c^ MS ), where 4 MS (p) -> 1.5p 2 = -23.2. In Fig. 2 r^(p = 27) is plotted as a 
function of a, c 2 near this saddle point, together with the contours of the constant r^ 3 '. 
We have calculated the NLO as well as the NNLO approximants of the quantity in the 
first square bracket of (0), with the unity subtracted, in the three chosen RSs and in the 
interval p G (18,28), which corresponds to the measured value of Mz = 91.4 GeV and 
in the interval e (50, 500) MeV. As the differences are tiny we normalize all our 
results to the NLO result in the conventional MS RS, and plot the relative difference, 



-W 



(RS) 



r '-rNLO (A _) L C 63 ) 

where i=NLO, NNLO and RS=PMS, ECH or MS. We draw the following conclusion 
from Fig. 3: 

1. The differences between PMS and ECH approaches are minuscule, about 0.1% at 
the NLO and totaly negligible at the NNLO, 

2. the difference between the PMS (or ECH) and MS approaches is 

• about 0.7% at the NLO and 

• about 0.3% at the NNLO (this comes from the ratio of the dotted and dash- 
dotted curves in Fig. 3). This documents the trend, observed in earlier works, 
that inclusion of the NNLO corrections diminishes the RS dependences and 
thus decreases the theoretical uncertainty. 
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Figure 3: The ratio r p , defined in the text. 
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3. the differences between the NLO and NNLO approximations amount to about 2% 
in the MS RS (as already estimated in Section 3) and to about 3% for PMS (or 
ECH). 

Given the current precision of the data (|5]), which translates into 8.5% accuracy on the 
couplant a(Mz), this implies that none of the mentioned differences is discernible. More- 
over, only the last difference, namely the effect of including the NNLO correction, has 
any chance of being seen in the data, but even this would require at least a factor three 
improvement in the precision of T/j. 



5 The dependence of Rz on and m t 

As already pointed out, massive quarks complicate the consistency conditions, discussed 
in Section 2, as their renormalized masses also 'run' and there is so far no generalization of 
the relations (|T^) available. We therefore estimate the effects of non-zero top and bottom 
quark masses within the MS RS. Recall that ([54]) are expansions in corresponding to 
five massless quarks, and thus all effects of finite mj, and m t are contained in the expansion 
coefliecients. 

First we deal with the effects of the finite x = Mz/m t in the first and third term 
of (|54|). In the first case the x-dependent terms of ri and r 2 vanish in the limit x — > 
and the leading, constant, ones correspond to rif = 5 in accordance with the decoupling 
theorem [[35|| . For x = 0.4225, corresponding to Mz = 91.187 and m t (Mz) = 140 GeV, as 
suggested in Ref. f2l|, R V ' NS + R A ' NS changes by the factor 3.7 10~ 5 , which is an effect 



two orders of magnitude smaller than the experimental error of Rz, estimated in (||), and 
one order of magnitude smaller than the effects discussed in the preceding section. It can 
therefore be safely neglected. 

The contribution of the last term in (|54"D, coming from the axial vector, flavour singlet 
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channel, where the familiar axial anomaly operates, is numerically nonnegligible. Despite 
the fact that the leading terms in t 2 and ts are x- independent, the nonvanishing of R A,S 
is due to non-zero m t , or, more precisely, to the effect of the difference m t — nib f\ This 
comes from the fact that axial couplings of quarks in weak dublets are opposite and, apart 
from mass effects, their contributions cancel. Consequently, for m t = both t% and £3, 
and thereby also R A,S , must vanish. The fact that this vanishing is not obvious in the 
small x expansions of ( fjif) is not surprising, as m t —>■ corresponds on the contrary to 
x — > 00. Note that the presence of powers of In a; in the coefficients t 2 and t 3 implies 
that a resummation of these logarithms is necessary |32| even in the limit x — > 0. As the 
contribution R A ' S can be interpreted as the top quark mass effect, we again determine its 
contribution in the MS RS only. For m t (M z ) = 140 GeV, R A > S = -0.00160. 

The effects of finite value of rrib/Mz have been studied in detail in [33]. Here we merely 
recall the form and numerical values of the leading contributions in the vector and axial 
vector, nonsinglet bb channels resulting from these effects: 



R 



V,NS 



R 



A,NS 



9v,b 1 + 



1 + 12 



m\{M z 



^(l-6^ + 0(a«; 



(64) 
(65) 



For m b (M z ) = 4.8 [0 and a (5) (M z ) = 0.037 the leading m b /M z contribution to R z 
coming from the axial channel thus amounts to 4.2 10~ 3 , whereas in the vector channel 
we get a mere 1.5 10~ 4 . 



6 The estimate of the still higher order terms 

As we saw in the preceding two sections, the NNLO corrections to are, compared to 
experimental accuracy of its measurement, tiny effects, and so it seems reasonably safe 
to stop at this order. On the other hand it is generally accepted that perturbative series 
in QCD do not converge, but represent merely asymptotic expansions to the full result. 
In such situation it is certainly useful to have at least some estimate of the magnitude of 
the so far uncalculated (and in the near future uncalculable) higher orders. 

An attempt in this direction has recently been made in P7fl , using the so-called 'im- 
provement formula' of Ref . |2j , which represents an approximation of the PMS optimiza- 
tion discussed in subsection 2.3. Its essence is to reexpand the PMS result optimized to 
the Nth order in powers of the couplant in any fixed RSQ and take the coefficient of this 
expansion at the (N+l)th order as an estimate of its true value. Instead of the PMS 
approach the ECH one of subsection 2.4 can be equally well used for this purpose. The 
resulting estimates are only slightly different. 

Here we outline the main steps of this method for the quantity fl3T)|), closely related to 
rV Consider first the Nth order partial sum of the perturbative expansion for a physical 



9 An alternative interpretation of this effect is discussed in Ref. |36[ . 

10 The optimized result is, of course, independent of the choice of this RS and is constructed from 
quantities up to the Nth order only. 
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quantity DP\ 

N-l 

= E d ia i+1 . (66) 

i=0 

Carrying out the optimization of ( |66|) according to either the PMS or ECH approaches 
leads to the optimized result, denoted below as (a op t)- If we now reexpand D^f (a opt ) 
in terms of the couplant a(RS) in a chosen RS we find 

DT(a opt ) = D N (a) + 5D&a N+1 , (67) 

where 

5DT = n N (di, a) - n N (d° p \ c° pt ) (68) 

give, according to ref. | 57] , the estimate of the coefficient in the chosen RS. For the 



three lowest orders the functions Qx(di, c») are given as |37 
fi 2 = d Q di(cx + di) 

Q 3 = d d 1 (c 2 --c 1 d 1 -2d 2 1 + 3d 2 ) (69) 

^4 = -^-(Sc 3 di + c 2 d 2 — 4:C 2 dl + 2cidid 2 — cxd 3 
+lAd\ - 2%d\d 2 + hd\ + 12did 3 ). 

These formulae can be derived from the following exact equations relating Qj to the 
coefficients dj,Cj and the RG invariants pf 



Pj c j _|_ 

J - 1 J - 1 d 



dj = -p^ - + ^. (70) 



Note that in order to evaluate Qj only di, q at lower orders % < j — 1 are actually needed! 

If the ECH approach is used for the optimization, the formulae ( |68| ) is particularly 
simple as d° pt = by definition and thus 

SDf™ = fi a (di,ci) 

<5D 3 ECH = tt 3 (4,4,ci,c 2 ) (71) 
<5.D ECH = n 4 (d 1 ,£/ 2 , d 3 ,cx,c 2 , c 3 ). 

The estimate ([71] ) is thus equivalent to the assumption that d N is dominated by the last 
term in (|70|) . Extensive discussion of this assumption and its consequences is given in 
Ref. |37| and is also related to Ref. [E3J]. 

Using the PMS approach the resulting estimate of d^ differs from ([FT]) by the pres- 
ence of the second term in (581), which does not vanish as in the ECH approach. How- 
ever, it was shown in ref. J37| that Q 2 (df MS , cf MS ) and fi 4 (<if MS , cf MS ) are small and 
fi 3 (<if , cf MS ) = 0. In the following numerical estimates only the ECH-based results are 
therefore presented. 

11 In this section we drop the specification of the number of quark flavours in the couplant a. 
12 The magnitude and therefore also the estimate of higher order coefficients depends on the choice of 
RS. We drop the argument 'RS' of a(RS) here. 
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n f 


r exact 


^estimate 
'2 


r est 


r estimate _ 


1 


-7.84 


-14.41 


-166 


-1750 


2 


-9.04 


-12.65 


-147 


-1161 


3 


-10.27 


-11.04 


-128 


-668 


4 


-11.52 


-9.59 


-112 


-263 


5 


-12.76 


-8.32 


-97 


67 


6 


-14.01 


-7.19 


-83 


330 



Table 1: The estimate of the so far uncalculated higher order coefficients for the quantity 
R(s) in the MS RS and using the ECH optimization procedure. 

There is one subtle point in the derivation of estimates for higher order coefficients 
of time-like quantities like ([1]) or (p6|). For instance, R(s) of ( flop is related to the so-called 
-D-function D(Q 2 ), defined primarily in the Euclidean region, via the dispersion relation 

B (« 2 »^ 2 f(jfw ds - (72) 

The knowledge of R(s) is in principle equivalent to that of D(Q 2 ), but as most of the 
optimization procedures or methods of higher order estimates do not commute with the 
functional on the r.h.s. of (0), we face the question as to which quantity to apply 
Eq.(71). This nontrivial problem is discussed in Ref. |37[| , the conclusion being that they 



should be applied to the quantities in the Euclidean region -, for instance, D{Q 2 ). Having 
obtained the estimates for higher order coefficients dj, the corresponding estimates for the 
coefficients rj of R(s) follow from the relations 





= di 






= d 2 - 


7T 2 b 2 


r-2 


12 






= d 3 - 


n 2 b 2 




4 




= G?4 — 


ix 2 b 2 


n 








4 



(2d 2 + -adt + -c\ + ca) + (73) 

The terms, proportional to powers of tt 2 , come from the analytical continuation of the 
couplant a(n) from the Euclidean region, where fi 2 < to the Minkowskean one, where 



fi 2 > 0. Taking into account the fact that in the MS RS we have |39 

di(MS) w 1.986-0.115/1/ 

d 2 (MS) w 18.244 -4.216^ + 0.086^ (74) 

77139 - 15099n f + 325n 2 f 

C2 < MS > = 9504 - 5 76, v ^ 

and using (FT]) we get the estimates, obtained originally in Ref. P7| , summarized in Table 
2 and valid for the MS RS.0 In order to get some feeling of the possible accuracy of these 

13 



Neglecting the terms of the light-by-light type, which violate the structure of (pq 
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estimates the above table also includes the results for the NNLO coefficient r 2 , for which 
the exact calculations are available. In the case of the coefficient only the estimate for 
the combination r 4 — r\c^ is presented, as the four loop /^-function coefficient C3(MS) is 
so far unknown. 

As the dominant contribution to ([!]) comes from the nonsinglet channel - first term of 
([54]) -, the above estimates are relevant for this quantity as well. For a(MS, Mz) = 0.037, 
rif — 5 and using the estimates of Table 1, we find that the terms r^a 4 and r^a 5 contribute 
approximately —3 10~ 4 and 8 10~ 6 respectively^. Note that while the latter contribution 
is entirely negligible, the former is of the same order as the RS uncertainty of the NNLO 
contribution. 

7 Summary and conclusions 

In the preceding sections we have analyzed various contributions to, and theoretical un- 
certainties of, the quantity (0). The results of these analyses are summarized in Table 2. 
All these numbers should be contrasted with the current experimental error of Rz, which 
is 5.9 10~ 3 . We see that there is a number of affects that are comparable to the current 
experimental accuracy of Rz, the most important of them being the NLO perturbative 
correction and, interestingly, the effects of finite b quark mass correction to the Born 
term in the axial channel. This, however, has nothing to do with QCD. On the other the 
hand, the data are not yet sufficiently precise to be sensitive to, for instance, the NNLO 
perturbative correction. Further improvement in the measurement of is clearly very 
desirable. 
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For the latter contribution the additional assumption, concerning C3(MS), had been made: C3 = c|/c. 
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type of contribution 


contributes to Rz 


1 


LO nonsinglet channel 


+62.2 icr 3 


2 


NLO nonsinglet channel 


+3.24 1(T 3 


3 


NNLO nonsinglet channel 


-1.08 10" 3 


4 


N 3 LO nonsinglet channel 


-5.1 10" 4 


5 


N 4 LO nonsinglet channel 


+1.3 10~ 5 


6 


m t in nonsinglet channel 


+3.7 10~ 5 


7 


m t in singlet channel 


-1.6 10" 3 


8 


smooth thresholds in a s 


+1.3 10" 3 


9 


mj) effects in Born term, axial channel 


-4.2 10" 3 


10 


rrib effects in LO term, vector channel 


+ 1.5 10^ 4 


11 


RS uncertainty at NLO 


-4.7 10" 4 


12 


RS uncertainty at NNLO 


-2.1 10~ 4 


13 


experimental error 


±5.9 10" 3 



Table 2: Summary of various contributions to, and uncertainties of, Rz- All numbers, 
except items 8,11,12 are correspond to MS RS. 
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